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ABSTRACT 

Anticipating new spectroscopic studies of quasar companions, we calculate expected 
emission line spectra for the companion galaxies. Under the assumption of isotropic pho- 
toionizing radiation, we predict equivalent widths for Ha, H/3, [O m] A5007, [O n] A3727, 
[Ne in] A3869, and Mg n A2798 lines. We consider 6 quasars, two of which (PKS 2135- 
14 and PKS 2349-014) have very close companions with projected separation < 4/i _1 
kpc. For these companions the calculated emission lines should be easily observable, 
with W\ ;> 20 A for the [O in], [O n], and Ha lines. Presence of strong emission lines 
from the companions occupying only a tiny solid angle would suggest that 1) quasar 
radiation is isotropic; 2) companion galaxies have measurable amount of gas, which can 
be assessed using the photoionization models; 3) host galaxies, if present, have been 
depleted of gas or do not interact with the companions. 



Subject headings: quasars: emission lines — quasars: general — quasars: individual 
(PKS 2135-14, PKS 2349-014, PHL 909, NAB 0205+02, HE 1029-1401, PKS 0405-123) 



1. Introduction 

There is increasing observational evidence that quasar activity is associated with compact 
groups of galaxies ( [Bahcall Chokshi 1991 ; Yee fc Ellingson 1993 , and references therein). Recent 



HST observations (Bahcall, Kirhakos, & Schneider 1995a, b, 1996; Bahcall et al. 1997| , hereafter 



BKSS) of nearby luminous (z < 0.3; My < —22.9) quasars revealed extremely close companion 
galaxies, within 2" from the QSOs (~ 4/i _1 kpc in projection). Table || summarizes the known 
parameters for the companions of six nearby quasars. The first three columns give quasar's name, 
its redshift and absolute magnitude. The other columns describe the quasar companions, including 
their redshifts, magnitudes, and projected separations from the quasar, d, both in arcseconds and 
kiloparsecs. The last column is the covering factor w, discussed in §|2[ 

The physical conditions in the interstellar medium of the companions are different from the 
Galaxy. Powerful ultraviolet radiation from the QSO, if isotropic, is able to ionize a large fraction 
of the gas. For a distant observer, the companion galaxies would look like giant HII regions. 
Therefore, spectroscopic observations of the companion galaxies of quasars directly test anisotropy 
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of the quasar ionizing radiation suggested by the implied AGN model (cf. for example Antonucci 



1993). The results of such observations may be important in understanding the quasar environment. 



i : i 

Stockton & MacKenty (1987) observed extended line emission around a sample of quasars in 
their narrow-band [O in] A5007 survey. In many cases, the emission regions are not confined to 
the nearest detected companion and instead show asymmetric structures, such as bridges and tidal 
tails. The Stockton-MacKenty results suggest that at least in some cases companion galaxies are 
experiencing violent interactions with the QSOs and that their gas may have been removed by 
interstellar shocks. However, very little observational evidence exists on the emission-line spectra 
of the close companions. Our work is intended to motivate future observations by showing that gas 
content and metallicity can be assessed by obtaining spectra of the companions. We calculate the 
emission line spectrum expected from a single gas cloud with various galactic parameters, which is 
illuminated by the hard UV radiation from a central QSO source. Comparing our models with the 
results of narrow band imaging or high-resolution spectroscopy yields estimates of the gas density 
and its distribution. 

We outline in §^ the calculations and the method that can be used to study the ISM in the 
companion galaxies. Next, we discuss in §|3]the results for individual companions. Finally, in §|], we 
speculate on what can be learned about quasars and their environment from the study of emission 
lines of the companion galaxies. 



2. Computing the emission lines 



We use the photoionization program CLOUDY , version 84.12 flFerland 1993 ) to calculate the 



emission line spectrum of a cloud of gas exposed under the quasar radiation. The density of the gas, 
njj, is assumed constant throughout the cloud. The shape of the incident QSO's spectrum is from 



Mathews fc Ferland (1987) and is normalized to the total quasar luminosity, Lqso- The underlying 
continuum of the companion is assumed to be the standard spectrum of a Sbc galaxy, as reported 
by |Coleman, Wu, fc Weedman (1980) , normalized to L comv . We explore five orders of magnitude 



in the gas density (values characteristic of a normal galaxy), and four orders of magnitude in the 
ionizing flux. The equivalent widths of the emission lines are evaluated relative to the outgoing 
continuum of the companion galaxy that includes the component of the incident radiation absorbed 
and re-radiated by the gas. For a similar exploration of clouds within, not outside, the AGN see 
Baldwin et al. (1995). 



The gas in the cloud will be partially (or fully) ionized by the quasar radiation and will resemble 
a giant HII region. The principle differences from galactic HII regions are the geometry of the cloud 
and the shape of the ionizing continuum. The companion galaxy usually occupies a small solid angle 
as seen from the quasar, and therefore only a fraction of the ionizing flux will hit the inner surface 
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of the cloud. This fraction is characterized by the covering factor 

where Q is the solid angle subtended by the cloud, S is the surface of the cloud closest to the 
quasar, and d is the distance to the quasar. 

Unlike a galactic HII region, the ionization region around the quasar does not have a definite 
boundary because of the power-law spectrum of the ionizing radiation flBahcall, Bahcall, &: Ko~ 



zlovsky 1971; Halpern & Steiner 1983). The electron density declines smoothly with radius and it 



is not possible to separate completely the ionized and neutral zones. Nevertheless, it is instructive 
to introduce a parameter analogous to the Stromgren radius of the HII region. Since CLOUDY 
imposes spherical geometry on the cloud, the ionized region is a spherical shell with surface area S 
and mean radial width ARs- 

a B n e n p bf j \ d J \ n H J 

where cub ~ 4 x 10~ 13 cm 3 s _1 is the hydrogen recombination coefficient, and n e and n p are 
the electron and proton densities, respectively. A filling factor, /, accounts for the non-uniform 
distribution of the gas in the cloud. Q(H) is the total rate of ionizing photons from the quasar. 
For the chosen shape of the quasar continuum, Q 55 = Q(H)/10 55 s" 1 = L Q50 /3.2 x 1O U L . The 
"Stromgren width", ARs, provides a useful estimate of the effective volume occupied by the ionized 
gas that determines the line emission strength from the cloud. 

Table | gives the values of ARs for our sample. Only two close companions have ionization 
zones comparable to their size, and in the rest of the galaxies the Stromgren zone is ionization 
limited (its size AR > ARs). The actual widths of the companions (column 4), as well as their 
covering factors zu, were estimated from the HST images. 

In the CLOUDY calculations we use solar abundances of elements. Our choice is supported 
by the recent result that the abundances of the heavy elements considered here (O, N, Ne, and 
Mg) differ by less than 50% in the solar neighborhood QEdvardsson et al. 1993| ), and by the current 
models of chemical evolution of galaxies (for example, [Pagel 1989 ) that predict slow evolution of 



metallicity over the last 5 — 10 Gyr. External galaxies, however, show a spread of metallicities. 
While including abundances as free parameters would significantly complicate our models, it is still 
necessary to investigate the extent to which varying metallicity of the gas may change the predicted 



spectrum. Interstellar dust can also modify the thermal balance in the clouds ( Shields fc Kennicutt 



1995[ ). We have conducted test calculations including the depletion of certain elements onto dust 



grains (as described in Baldwin et al. 1991| ) and varying metallicity parameter in a wide range of 
O.1Z < Z < 5Z Q . 

At very low metallicity (Z ~ O.IZq), the forbidden line emission (for example, [O in] ) increases 
as Z rises, while the hydrogen Balmer line intensity is essentially constant. Increasing metal 



- 4 - 



abundance leads to more efficient cooling and correspondingly lower electron temperature. At 
about solar metallicity the temperature effect starts to dominate and both [O III] and H/3 decline 



even as Z increases (McCall, Rybski, & Shields 1985). Including dust grains in the calculations 



gives rise to the opposite phenomenon. The main effect of the grains is depletion of metals from 
the interstellar medium ( Shields fc Kennicutt 1995| ) . This reduces cooling efficiency of the gas and 



leads to higher temperature than in the dust-free case. Therefore the variation of the predicted 
emission lines with Z is reduced. 

The test calculations have been performed for all quasar companions considered in §||. In the 
range 0.3 < Z/Zq < 3, the intensity of the [O in] line varied by an average of 45% for the two 
close companions and by 53% for the other companions with lower ionization level. In the range 
0.1 < Z/Zq < 5 for the close companions, [O in] was reduced by about 90% compared to its 
value at Z ~ Zq. The [O III] /H/3 line ratio changes more slowly than the [O in] intensity at high 
metallicity because of similar drop in H/3. In all calculations involving the two close companions the 
average electron temperature was in the range T e = 8 x 10 3 — 3 x 10 4 K, and for the low-ionization 
companions T e = 6 x 10 3 — 2 x 10 4 K. Prom these tests we conclude that choosing solar abundances 
maximizes the forbidden line emission from the quasar companions. 



2.1. Distribution and density of gas in the close companions 

When companions are very compact or very close to the quasar, such that AR < ARs, 
essentially all gas is ionized and the intensity of the emission lines is directly proportional to the 
effective volume occupied by the gas, V = S x AR x /. The unknown parameters are the gas 
density n# and the filling factor /. We have computed a two-dimensional set of ionization models, 
varying n# from 10 -2 to 10 3 cm -3 , and / from 1 to 10 . 

Figure |l] shows the contours of the equivalent widths for the Ha, H/3, [O in] A5007, [O n] A3727, 
[Ne in] A3869, and Mg n A2798 lines, corresponding to the rest frame of the cloud. The [O in] line 
illustrates the transition from the matter-bounded to ionization-bounded models. On the left side 
of the plot ARs > AR and the contours of W\ coincide with the lines of constant Stromgren radius 
(/ oc n^ 2 ). From the simple three- level model of an oxygen atom we can predict Wa([0 hi]) cx n^jV. 
For large densities, ARs drops below AR and the line intensity becomes essentially independent 
of the amount of gas, /. Also as density increases, the ionization level of oxygen decreases and the 
[O n] line becomes stronger. 

In the Ha, H/3, and [Ne in] panels there is a plateau on the high-density side where the 
equivalent widths are essentially constant. Here all of the ultraviolet photons are trapped inside 
the cloud, and the ionization zone is limited by the quasar intensity. For example, W\(Hf3) cannot 
rise much more than 30 A for any combination of the gas density and the filling factor. 



Note that [O in] A4959 line is also present if [O in] A5007 is observed. We do not include this 
line because of the approximate relation Wa(5007) ~ 3 x W\(4959) for T ~ 10 4 K (Osterbrock 
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1989). If the two lines are not resolved, the equivalent width of the [O m] A5007 line should be 
corrected accordingly. 

The intersection of the isocontours of all the observed W\ gives the point on the f — nu diagram 
that corresponds to the ionization model best reproducing the observations. Besides an estimate 
of tih, we obtain also the full emission-line spectrum of the object (within the simple, constant 
density model). Therefore we predict the intensity and equivalent widths for other strong lines, and 
these in turn can be tested by future observations. Given good data from several lines (including 
upper limits), it seems possible to construct a model with a realistic gas distribution covering a 
range of densities nn- 



2.2. Density of gas in the ionization-bounded companions 

If the companions are sufficiently far from the quasar, the strength of the emission lines is 
limited by the number of ionizations. At high enough density, the effective volume of the ionized 
region is V = S x AR$ X / and the equivalent widths of the hydrogen lines are proportional to a 
single ionization parameter 

- _ Lqso 



■^comp 



(3) 



This parameter is, in principle, directly observable. Values of 5 for the companion galaxies in our 
sample are given in the last column of Table |2[ It could be related to the usual ionization parameter 
U: 3 = U X n/f(cm" - 3 )5 , (kpc 2 )(L comp /8.8 x 10 9 Lq)- 1 . The two paramet ers are similar in value 
but H has the advantage of being independent of the gas density. This allows us to separate the 
intensity of the incident flux from the parameters intrinsic to the companion. 

Figure ^ shows the results of our calculations for the same emission lines, now as function of 3 
and riH- The two hydrogen lines, Ha and H/3, demonstrate the simple variation of the equivalent 
widths. When the density is large, W^(Ha, H/3) oc 3 almost independently of n#. On the other 
hand, for nu < 1 cm~ 3 the ionization is complete and W\(Ha,}i(3) oc n H . The corners of the 
contours lie on a single line Ai? = AR$ (i.e. 3 oc n H ). Note that both permitted and forbidden 
lines show a general break along that line. 

Behavior of the oxygen lines can be understood using the three-level model of the oxygen 
atom. The dominant excitation mechanism is collisions with electrons unless the radiation field is 
so strong that recombination to excited levels becomes important. For the [O in] line, collisional 
excitation prevails if n(0 +2 )/n(0 +3 ) > 0.0025, which is always satisfied in our models. We have 
also checked that the quasar radiation is strong enough that photoionization of neutral oxygen 
atoms wins over the charge exchange with protons (see Spitzer 1978; ) if n p <^ 76 Q55 (d/5 kpc)~ 2 
cm -3 . 



For low densities ionization is so high that the ion abundances of the heavy elements are 
controlled by recombination from the next ionization state. As a result, all equivalent widths 
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follow the approximate relation W\(0, Ne, Mg) oc n^H -1 . In the limit of a weak radiation field (or 
high density), the ion concentrations scale differently and so do the equivalent widths of the lines. 
The ionization zone is limited by Ai? s , and Wa([0 ii]) oc E 2 nJ I 1 , W\([0 m]) oc H 3 n^ 2 . 

Figure [2] allows us to find the best ionization model describing a set of the observed emission 
lines. The ionization parameter of the cloud could be determined directly from the observations 
(see eq. [||). Ideally the contours of the equivalent widths should all intersect in one point with 
the horizontal line of the known 5, and even one measured line would suffice to obtain the value 
of the gas density. However given the uncertainty of measurements, it is desirable to have several 
lines in order to better constrain the estimate of n#. 

How does the observational uncertainty in W\ propagate to the inferred values of n#? We 
have considered a range of values of W\ and obtained the corresponding range of allowed densities. 
Our tests show that the relative error in log nu is of the same order as that in log W\ ( [O in] , [O n] ) , 
with Alogn# ~ 1.3 A log W\ for a wide region on the 5 — nu diagram. Therefore, observational 
errors are not strongly enhanced in the estimate of the gas density. 



3. Discussion of individual galaxy companions 

In this section, we apply the results of our calculations to the individual companions listed in 
Table [l]. In calculation of luminosities of quasars and their companions we use a solar V magnitude 
of My iQ = 4.84 QWorthey 1994] ) and the Hubble constant H = lOO/i km s _1 Mpc" 1 . Note that 



choice of Hq does not affect the predicted equivalent widths because the ionizing fluxes from the 
quasars are invariant to distance scaling. 

To account for the third spatial dimension, we correct the projected separations of the com- 
panions from the quasars by a factor of y/3/2 ~ 1.22. These larger values are used in calculation 
of the ionization models. The second column of Table ||| gives the corrected separations for all of 
the companions. 



3.1. Close companions 

3.1.1. PKS 2135-14 

I .=1 

Stockton (1982) presented a spectrum and a detailed discussion of the closest companion to 
PKS 2135-14. Its compact nature is best revealed on the HST image (BKSS), where the companion 
is unambiguously distinguished from the nebulosity surrounding the quasar. Based on the spec- 



trum of the companion, Canalizo & Stockton (1997) showed that it is actually a foreground star. 



Therefore, the narrow line emission is produced by the gas in the fuzzy nebulosity extending over 
more than 8/i _1 kpc. Our ionization models describe the state of this extended gas. 
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We read the equivalent widths of the [0 III], [O n], and H/3 lines by visual inspection of the 
spectrum in Stockton (1982) . Table [3| shows our estimates of W\\ uncertainty in the values should be 



taken to be at least 20%. We construct the best fitting photoionization model using the calculations 
described in § |2.1| . Figure [3| shows contours of W\ for the three measured lines. One striking feature 
of the plot is that all three contours run reasonably close along the line corresponding to the 
constant Stromgren radius (/ oc n 2 H ). Therefore any model along the line is, in principle, capable 
of reproducing the observed emission. However, the isocontours never cross in one point, suggesting 
that more than one constant density model is required for the adequate fit of observations, provided 
the estimated equivalent widths are correct. 

We considered a simple combination of high and low density clouds (see Table ||). The high 
density clouds are needed to reproduce the [O n] emission, since most of the oxygen is doubly 
ionized at small densities (jih < 3 cm -3 ) at the apparently small observed separation from the 
quasar. The low density clouds are likely to be present in the large extended nebulosity. We stress 
that the choice of high and low density clouds is not uniquely set by the available data. 

The two models together produce the right amount of the [O n] emission, 23% less of the 
[O in], and 50% more of the H/3 than is estimated from the observed spectrum. The Ha line should 
be conspicuous, with W\ ~ 24 A, but it lies outside the spectral range of Stockton's data. The 
[Ne v] A3426 line becomes moderately strong if we take any lower density for the model a. However 
it is barely detectable in the spectrum, and the estimated upper limit on its equivalent width is 
about 3 A. Thus the gas density of the model a is limited to nn ;> 2.5 cm -3 . On the other hand, 
a larger filling factor (/ > 10~ 2 ) leads to a significant overestimate of the H/3 line while still not 
producing enough of the [O in]. 



Another companion is located 6" to the East from PKS 2135-14. Canalizo fe Stockton (1997) 
confirmed that it is a galaxy at the redshift of the quasar, however its spectrum does not show 
strong emission lines. 



3.1.2. PKS 2349-014 

The closest companion of PKS 2349-014 exhibits a tidal interaction with the host QSO (see 
Bahcall et al. 1995a). Its redshift has been confirmed by Miller, Tran & Sheinis (1996). The 



estimated three-dimensional distance from the quasar is 4.2/i _1 kpc. There are also apparent tidal 
tails and a large nebulosity surrounding, but offset from, the quasar. If the quasar has a host galaxy 
similar to the Milky Way, the companion would inspiral due to dynamical friction in 7 x 10 8 yr. 

We have computed two illustrative models with plausible gas densities of nn = 1 cm~ 3 and 
4 cm -3 , respectively. We expect that the highly ionized, and therefore heated, gas would tend to 
be evenly distributed over the volume of the companion. Because of the faint continuum of the 
companion, the predicted equivalent widths in the second model are very large (see Table ||). For 
example, W^([0 ill] ) ~ 250 A. Spectroscopic observations of the companion galaxy are necessary 
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to constrain the ionization model and to estimate the present amount of gas. 

Because of the extreme proximity of the companion, radiation pressure from the quasar may 
affect the distribution of the gas. We use CLOUDY calculations to estimate the radiative accel- 
eration on a hydrogen atom in the companion. We assume an isothermal distribution of mass of 
the companion and adopt the mass-to-light ratio of M/Ly = 1 characteristic of moderately young 
stellar population. We find that for the models I and II, gravity of the companion dominates the 
radiation pressure out to 18 kpc and 8 kpc, respectively, much farther than the actual extent of 
this object. Therefore, the gas clouds are not driven out of the companion. We have checked 
that the momentum transfer due to absorption of ionizing radiation by hydrogen atoms (for exam- 
ple, [Haehnert 1995 ) is at least an order of magnitude weaker than the force due to resonant Lya 



scattering. 



3.2. Ionization-bounded companions 

3.2.1. NAB 0205+02 

NAB 0205+02 from the survey by Bahcall, Bahcall, & Schmidt (1973) has a mysterious com- 



panion at a projected separation of 22h kpc; this companion is very compact but is bright in 



the [O in] line. Stockton & MacKenty (1987) estimated the equivalent width of [O in] line of at 
least 400 A from the narrow line image. On the other hand, the object was not detected at all on 
their line-free continuum image. Recent HST F606W broad band observation of the companion 
(BKSS), which includes the [O m] line, shows that the image consists of two tiny clumps connected 
to each other by a bridge. There is no evidence on the HST image that the quasar host extends all 
the way to the companion, suggesting that the gas in the object is not affected by direct interaction 
with the host. The line emission in the companion may, in principle, be caused by photoionization 
or collisional excitation by sources internal to the companion. We assume here that it is not the 
case and discuss instead the situation when photoionization by the quasar is the dominant source 
of energy for the emission lines. 

The Stromgren radius for the companion (see Table ^) is less than the apparent extent of 
the object (~ 3.5 kpc in NS direction, and 0.5 kpc in SW direction), so that the gas will be only 
partially ionized for nu > 0.3 cm -3 . The amount of [O in] emission is therefore determined by the 
ionization parameter S, which in turn depends on the covering factor. Assuming that the size of 
the object along the line-of-sight is the same as in the NS direction, we get an upper limit on w 
of 1.4 x 10 -3 , or S ~ 1. Using the estimated three-dimensional separation from the quasar of 26.5 
kpc, we found that no simple ionization model, with the ionizing source centered on the quasar, 
yields an equivalent width of the [O in] line as large as is observed. The smallest covering factor 
needed to produce Wa([0 III]) ^ 400A is w = 2.4 x 10~ 3 , and the gas density in this case is 0.3 
cm -3 . 
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Can the emitting gas be distributed more broadly than the HST image shows? Apparently, 
the narrow-line image of Stockton & MacKenty (1987; their Figure 2e) is large but this could be 
the effect of lower resolution relative to HST . Nonetheless, we present the ionization model with 
w = 2.4 x 10~ 3 that best fits the observations. If the model is confirmed by new spectroscopic 
data, a better understanding of the distribution of gas and stars in the companion will be required. 
If not, a much more elaborate ionization model will be needed to describe this unusual object. 

The one-cloud ionization model predicts large equivalent widths for virtually all of the emission 
lines listed in Table f|. The [O n], Ha, H/3, [Ne in], and Mg n lines have calculated W\ in excess of 
40 A. Also included in Table ^ are strong [N n] and ITy lines. Our simplest model gives an upper 
limit (for / = 1) on the total mass of the gas in the companion of M gas ^8 x 10 7 Mq. 

The faint continuum of the companion suggests the possibility that most of the light might be 
produced by thermal bremsstrahlung rather than starlight. The bremsstrahlung luminosity of the 
companion would be ( Psterbrock"T989|) 

L ff = 1.42 x 10' 27 g ff T 1/2 n e niZ 2 V erg s" 1 ps 1.5 x 1O 7 L T 4 1/2 n e niZ 2 , (4) 

where gjf ~ 1.3 is the free-free Gaunt factor, rij is the density of ions of charge Z, and V ~ 10 
kpc 3 (for w = 2.4 x 10~ 3 ) is the volume of the gas. The observed companion luminosity is 
^cornp = 1-3 x 10 8 Lq (BKSS), and therefore one needs a density of n e > 3 cm -3 to produce the 
right amount of emission. This in turn requires a larger covering factor of w > 6 x 10~ 3 in order 
to match the observed W>,([0 III]) with the gas density larger than 3 cm -3 . As noted before, the 
covering factor can hardly be larger than (2 — 3) x 10~ 3 . Also, self-gravity of the gas will not be 
strong enough to resist the radiation pressure force, as indicated by the CLOUDY calculations. 
Therefore the continuum of the companion cannot be accounted for by thermal bremsstrahlung 
only. Most likely, the detected emission arises from a combination of emission lines and faint 
starlight. 



3.2.2. PEL 909 



Galaxy F in the Bahcall, Kirhakos, &c Schneider (1996) list of objects around PHL 909 is 65ft,~ 1 
kpc away from the quasar. [Boroson (1996) has estimated the equivalent widths of the [O n] and 
[O Hi] lines as 10 A and 2 A, respectively. 

Figure shows the corresponding contours of the lines as function of H and n#, using the 



calculations from §[2^. The intersection of the two lines gives a slightly larger covering factor than 
is estimated from the HST image, so that we had to lower the gas density in order to produce 
the observed W\. The best ionization model is obtained with uh = 0.8 cm~ 3 for w = 5 x 10 -3 . 
Because of the large distance to the quasar, the ionization flux is not very powerful and there are 
no strong lines predicted. Besides the oxygen lines, there could be detectable Ha and [N n] A6584 
lines both with W\ ~ 4 A. 
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3.2.3. HE 1029-1401 

1 ; — 1 

Wisotzki (1994) showed that the quasar HE 1029-1401 is surrounded by a loose group of 
galaxies with strong emission lines. Table |] gives the equivalent widths for the [O III], [O n], 
and H/3 lines for the closest galaxy with a redshift similar to the quasar. The companion is at 
deprojected distance of about 150/i _1 kpc, and the covering factor is small (w ~ 6 x 10~ 4 ). Even 
with w = 10~ 3 and nu = 1 cm -3 , photoionization by the quasar is too weak by an order of 
magnitude to produce the measured lines (Table Q). 

The observed emission can be explained by star formation within the companion galaxy. The 
continuum level at H/3 (A = 486lA) is related to the total luminosity L COUiV using the standard 
spectrum of an Sbc galaxy from Coleman et al. (1980); we find XL\(K(3) = 0.93L comp . The 
observed equivalent width of the H/3 line gives then the intensity of the line, L(H/3) = 2.1 x 10 40 /i~ 2 
erg s _1 . Using Kennicutt's (1983) relation for the star formation rate (SFR) of spiral galaxies to 



the luminosity of the Ha line and assuming L(Ha) ~ 3 x L(H/3), we find SFR ~ 0.6 Mq yr 
which can be achieved by a normal size spiral galaxy. 



3.2.4. PKS 0405-123 



The companion galaxy to PKS 0405-123 has been observed by |Marr fc Spinrad (1985) , and its 



redshift has been confirmed by Ellingson fc Yee (1994) . The spectrum of the galaxy has a strong 
[O 11] emission line, with the estimated equivalent width of IUx([0 11]) ~ 12 A. The luminosity 
of the companion (My = —20.2) suggests that it may be a large spiral galaxy. There is no HST 



observation of this quasar field and the available image from Marr fc Spinrad (1985) does not permit 
robust determination of the covering factor of the companion. We estimate the size of the galaxy 
to be 10 kpc, and assume it is almost face-on towards the quasar. This gives w ~ 3 x 10" 3 . Using 
Figure |2j, we estimate the gas density to be about 8 cm -3 . The ionization model predicts some 
measurable emission in Ha, Mg 11 and [N 11], with W\ ;> 5 A. If the companion is inclined with 
respect to the quasar, w will be smaller and the ionizing flux from the quasar will not be strong 
enough to produce the observed [O 11] emission. 



4. Summary 

The primary goal of this paper is to encourage new spectroscopic observations of quasar com- 
panions in order to learn more about the companion galaxies, isotropy of the quasar emission, and 
the relation between the galaxies and the quasars. We have calculated the equivalent widths of 
the strongest emission lines expected to be produced by quasar photoionization of close galactic 
companions and applied our models to the six companions with the observed lines. 

In our exploratory calculations we have assumed constant density gas clouds and solar abun- 



- 11 - 



dancies of the elements. Choice of the solar metallicity maximizes the forbidden line emission. For 
the close companions of PKS 2135-14 and PKS 2349-014, the predicted equivalent widths are very 
strong, W\ ^ 20 A. Our results suggest (see Figures Q and |, and Tables | and |) that a narrow 
range of gas densities of nn around 1 cm~ 3 may produce the observed emission lines. In three 
cases (the two close companions mentioned above and the exotic companion of NAB 0205+02) the 
models constrain the amount of gas to be of order M gas ~ 10 7 — 10 8 Mq. 

Determining the amount of gas in or around the close companions is important for understand- 
ing the dynamical interactions of quasars with their environment. HST images (e.g., Bahcall et al. 
1995a and BKSS) show clearly that, at least in some cases, quasars are interacting violently with 
their galactic companions. The long tidal tail observed near PKS 2349-014 suggests the presence 
of a massive host galaxy. The small quasar-companion separation of ~ 4/i _1 kpc indicates that the 
companion is well inside the typical scale of a large galaxy. Gas shocks in the interstellar medium of 
the host might be expected to remove all of the gas from the small companion. Since the radiation 
pressure of the quasar is not strong enough to push away all the gas of the companion, absence 
of the strong emission lines predicted by our isotropic photo-ionization models could suggest an 
anisotropic ionizing flux or shock interactions with the massive host.[] On the other hand, if strong 
photoionized emission lines are observed, a possible explanation might be an almost "naked" central 
quasar, which accretes gas from disrupted dwarf companions. This possibility has been suggested 
recently by Bahcall et al. (1995a) and by Fukugita & Turner (1996). 

The degree of anisotropy of the quasar ionizing radiation, which is important for understanding 
the unified AGN model (e.g. Antonucci 1993| ), can be tested directly. If the continuum emission 



from the quasar is concentrated in narrow beams, the small solid angle of the companion galaxies 
makes them unlikely to be ionized. The companion of NAB 0205+02 presents a particularly strong 
test. Its tiny covering factor (w ~ 10~ 3 ) guarantees that the companion will not be ionized unless 
the quasar radiation is isotropic. Spectroscopy of that object can determine if the line emission is 
consistent with photoionization by the quasar. 

The gas in the host galaxy itself should be strongly ionized by the quasar and all of the emission 
lines considered in this paper are expected to be observable. For the simplest model of a disk galaxy 
with the constant gas density, the effective radius of the ionization zone is Ri on ~ 4.0 (Ln/Jn^) 1 ^ 



kpc, where L\\ = Lqso/W i Lq. As noted by Begelman (1985) and Chang, Schiano, fe Wolfe 



(1987) , total depletion of gas and dust occurs in the inner few kiloparsecs. Thermal expansion of 
the gas in the direction perpendicular to the galactic plane leads to an amorphous envelope of hot 
gas around the quasar extending up to several tens of kpc. [Stockton fc MacKenty (1987) found 



such large emission envelopes to be ubiquitous around quasars. Therefore the period of quasar 



An alternative explanation for the lack of line emission in the companions (proposed by J. Goodman) could be 
a very short lifetime of the quasar. If the central engine is powered by the disrupted material of a single star and 
lasts for only several years, the ionizing radiation from the quasar may have not yet reached the companions some 
10,000 light years away. 
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activity should be longer than the timescale for thermal expansion and subsequent cooling of the 
gas, or ~ 10 8 years. 

Close companions to three other quasars, PKS 1302-102, 3C 323.1 (=Q1545+210) and PG 
1202+281, have been detected with HST ( [Bahcall, Kirhakos, fc Schneider 1995b| ) but no spectra 
are currently available. [Stockton fc MacKenty (1987) found an extremely strong [O in] emission 



around Q1545+210, but they could not unambiguously associate the emission with the companion 
7 kpc away from the quasar. We would like to attract reader's attention to the 6 companions 
discussed in this paper and the other 3 objects mentioned above. They present a new interesting 
opportunity to understand quasars. 
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and [O in] emission lines on his spectrum of galaxy F in the field of PHL 909. We thank J. Bahcall, 
B. Draine, G. Ferland, J. Goodman, J. Ostriker, M. Rees, and L. Spitzer for helpful discussions. 
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Table 1. Photometric Parameters of Companion Galaxies' 



Quasar Companion 





z 


My 


z 


My 


d(") 


d{h~ l kpc) 


w b 


PKS 2135-14 


0.200 


-23.5 


0.200 (?) 


-19.4 


2.0 


4.2 


£ 2 x 10~ 2 


PKS 2349-014 


0.173 


-23.4 


0.173 


-17.7 


1.8 


3.4 


7 x 10~ 3 


NAB 0205+02 


0.155 


-23.0 


0.155 


-15.4 


12.0 


21.6 


£ 1.4 x 10~ 3 


PHL 909 


0.171 


-22.9 


0.169 


-20.1 


28.2 


53 


5 x 10~ 3 


HE 1029-1401 


0.086 


-23.2 


0.086 


-18.5 


112 


120 


6 x 10~ 4 


PKS 0405-123 


0.574 


-24.3 


0.570 


-20.2 


13 


49 


3 x 10~ 3 



a Using H = 100/t km s _1 Mpc -1 , q = 0.5. 

b Covering factor w of the companions, estimated from their HST images. 



Table 2. Ionization Parameters of Companion Galaxies 





df (h" 1 kpc) a 


AR S (kpc) b 


AR (h- 1 kpc) 




PKS 2135-14 


5.1 


1.5 


3 


0.9 


PKS 2349-014 


4.2 


2.3 


1 


1.3 


NAB 0205+02 


26.5 


4.1 x 10~ 2 


0.5 


1 


PHL 909 


65 


5.9 x 10~ 3 


9 


0.07 


HE 1029-1401 


147 


1.6 x 10" 3 


6 


0.05 


PKS 0405-123 


60 


2.7 x 10~ 2 


10 


0.1 



Reparation from the quasar corrected for the projection effect; d! = d x ^/3/2. 

b For Tin = 1 cm -3 , / = 1. The Stromgren radius scales as AR$ oc nj^f^ 1 (see eq. |2[). 
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Table 3. Emission lines for close companions 



PKS 2135-14 PKS 2349-014 







Model a 


Model b 








W X (A) 


[1] 


(low-density) 


(high-density) 


total 


Model I 


Model II 


[0 in] A5007 


69 


9.8 


43.3 


53 


53 


247 


H/3 


5 


2.0 


5.6 


8 


8 


24 


[0 n] A3727 


8 


< 1 


7.9 


8 


< 1 


19 


Ha 




6.3 


17.4 


24 


26 


73 


Mg II A2798 




< 1 


3.4 


3 


< 1 


25 


[Ne in] A3869 




< 1 


3.8 


4 


4 


22 


[Ne v] A3426 




2.3 


< 1 


2 


25 


6 






Parameters of ionization 


models 






tlh (cm -3 ) 




3 


45 




1 


4 


/ 




10~ 2 


io- 4 




1 


1 


M gas (M Q ) 




> 5.4 x 10 6 


8.1 x 10 5 


6.2 x 10 6 


2.4 x 10 7 


9.6 x 10 7 



References. — [1] [Stockton (1982] . 
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Table 4. Emission lines for ionization-bounded companions 





NAB 0205 


PHL 909 


HE 1029 




PKS 0405 


W A (A) 


[1] 


Model 


[2] 


Model 


[3] 


Model 


[4] 


Model 


[0 in] A5007 


^400 


415 


2 


2 


39 


< 1 




< 1 


[0 n] A3727 




118 


10 


9 


52 


6 


12 


12 


H/3 




48 




1 


13 


1.5 




2 


Ha 




147 




4 




5 




6 


Mg ii A2798 




94 




3 




3 




5 


[Ne in] A3869 




42 




< 1 




< 1 




< 1 


[N ii] A6584 




32 




4 




4 




6 


H 7 




27 




< 1 




< 1 




1 






Parameters of ionization 


models 








nji (cm -3 ) 




0.3 




0.83 




1 




7.5 


M gas (Mq) 




< 8 x 10 7 




< 10 10 








< 4 x 10 10 



References. - [1] |Stockton fc MacKenty [1557% [2] |Boroson (1996)1 ; [3] |Wisotzki (1993% [4] 
Marr fe Spinrad (1985] . 
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10" 2 0.1 1 10 10 2 10 3 0.1 1 10 10 2 10 3 
n H (cm -3 ) n H (cm -3 ) 



Fig. 1— Contour plot of equivalent widths for the Ha, H/3, [0 in] A5007, [O n] A3727, 
[Ne in] A3869, and Mg n A2798 emission lines as function of the filling factor, /, and density, 
nn, of gas in a close companion. Solid lines are the contours of W\ = 1, 10, 100 A, with the 
numbers labeling log W\. Dotted lines represent 0.25 dex steps. The other fixed parameters of the 
calculation are: Lqso = 2 x 10 11 Lq, L com p = 5 x 10 9 Lq, d = 3.2 kpc, AR = 2 kpc, w = 0.02. 
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10" 2 0.1 1 10 10 2 10 3 0.1 1 10 10 2 10 3 
n H (cm -3 ) n H (cm -3 ) 

Fig. 2.— Contour plot of equivalent widths for the Ha, H/3, [0 III] A5007, [O n] A3727, 
[Ne in] A3869, and Mg n A2798 emission lines as function of the ionization parameter, H, (eq. 
[||]) and density of the gas, njj, for the ionization-bounded quasar companions. Solid lines are 
the contours of W\ = 1, 10, 100, and 1000 A, with the numbers labeling logW>. The dots rep- 
resent 0.25 dex steps. The other fixed parameters of the calculation are: Lqso = 2 x 10 11 Lq, 
L comp = 2 x 10 8 Lq, AR = 1 kpc, / = 1. 
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10~ 2 0.1 1 10 10 2 10 3 

n H (cm" 1 ) 

Fig. 3. — Example of estimating the filling factor / and density Uh of gas around PKS 2135- 
14, using the map of ionization models from Figure |]. Solid curves are the isocontours of the 
observed equivalent widths for the three emission lines: Wa([0 hi]) = 69 A, W\([0 11]) = 8 A, and 
W\(H(3) = 5 A. Intersection of all the lines gives the best ionization model. In this case, however, 
[O ill] and H/3 contours are almost parallel to each other with some offset from [O 11]. 
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Fig. 4. — Isocontours of the emission lines for the companion of PHL 909: W>([0 11]) = 10 A and 
Wa([0 hi]) = 2 A. Intersection of the two lines gives a higher density than in the best-fitting model 
(njj = 0.83 cm -3 ) because of the smaller observed 3 parameter. The dashes outline the region in 
this two-dimensional space that is allowed by uncertainty in W\ of 30%. 



